Preprint typeset in JHEP style - HYPER VERSION 



Y-system and Quasi-Classical Strings 



Nikolay Gromov 



DESY Theory, Hamburg, Germany & II. Institut fur Theoretische Physik Universitdt, 
Hamburg, Germany & 
St. Petersburg INP, St. Petersburg, Russia 
E-mail: tiikgromov@gmail.com 



RE 



Abstract: Recently Kazakov, Vieira and the author conjectured the y-system set of 
equations describing the planar spectrum of AdS/CFT. In this paper we solve the Y- 
system equations in the strong coupling scaling limit. We show that the quasi-classical 
spectrum of string moving inside AdS^ x S 1 matches precisely with the prediction of the 
^-system. Thus the y-system, unlike the asymptotic Bethe ansatz, describes correctly the 
spectrum of one-loop string energies including all exponential finite size corrections. This 
gives a very non-trivial further support in favor of the conjecture. 
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1. Introduction 

Since the discovery of AdS/CFT correspondence [l], |2| there has been a significant progress 
in solving the maximally super-symmetric Yang- Mills theory in four dimensions, mainly 
due to the integrability found on both sides of the duality || |], |5[ |7j. The problem of 
finding the anomalous dimensions of all local operators in planar limit was thus brought 
very close to its complete solution. Recently the anomalous dimension of the simplest 
operator was computed numerically for a wide range of 't Hooft coupling A |J. 

The approach used there is based on the Y-system for AdS/CFT conjectured in B. 
The main ingredients leading to the conjecture are the asymptotic Bethe ansatz equations 
(ABA) |l(], pj]] , describing the spectrum of the anomalous dimensions of local operators 
with infinitely many constituent fields, the concept of the "mirror" double Wick rotated 
theory invented in [^] and then explored in details in [13], Liischer formula for the finite 
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Figure 1: T-shaped "fat hook". It defines the interactions between Y-functions in the AdS/CFT 
Y-system equations. 



size corrections |E4| and our experience with relativistic theories, where similar equations 
are constantly appearing in the thermodynamic Bethe ansatz (TBA) approach [15, p]]. 



This conjecture passes some nontrivial tests - in Q the 4-loop perturbative result [17] 



was reproduced. More recently a comparison was also made at 5-loops in [18]. In [^, ^(j] it 
was also shown to be consistent with the TBA approach and finally its numerical solution || 
indicates that the strong coupling asymptotic agrees with the string prediction |^] for the 
simplest Konishi operator 1 . Apparently more tests of this conjecture are indispensable. 
In this paper we compare the results of the quasi-classical string quantization with the 
prediction of y-system thus providing the first analytical test of the conjecture to all 
orders in finite size wrapping corrections. As we will see below this check involves numerous 
remarkable identities and miraculous simplifications and probes most of the Y-system in 
great detail. Therefore, the match we will observe leaves little room for doubt about the 
validity of the Y-system, at least in the strong coupling limit. 

The Y-system of |J is an infinite system of simple functional equations 2 

Y (*+±\Y ( Z M- (l + *W*))(l + ^, S -i(*)) n n 

as[ + 4e)W + l/ Ya+lyS (z))(l + l/n-U*)) ' 1 j 

The indices a, s belong to a T-shaped lattice (see Fig.|l|). One should replace Y as with indices 
outside the lattice by or oo so that they disappear from the equations. Note that the 
Y-system equations only involve the pair of indices a, s and its immediate neighbors only, 
a rb 1, s and a, s rb 1; in this sense these equations are described as "local". In accordance 
with Fig.g we also use the following notations 

{Ym a , Y/\ a , Yq s , Y®, Y®} = {Y a , , Y aA , Y l>s , Y hU Y 2j2 } . (1.2) 



lr rhe sub-sub-leading coefficient does not agree with the quasi-classical string results of pi| . However, 
the Konishi state is one of the lowest states and the quasi-classical quantization may not be applicable for 
this case. 

In this paper we use a rescaled rapidities z = where g — more convenient for our purposes. 
3 In the original derivation of the Y-system equations the Y-functions are associated with densities of 
several bound states of the mirror theory The symbols used as subscripts in this new notation were 

introduced in [[19| to indicate which type of bound-state they originate from. 
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To have a unique solution this system should be supplemented with an additional 
"non-local" equation [19] 1 



oo r?(+) r(-) 

log = £ (n^ - b^A * io g (i + y.j + bg ^ , (i.3) 



a=l 



and particular boundary conditions at infinite a and s. Here, 7Z and B are some ker- 
nels defined in Appendix A and * stands for convolution. This additional data encodes 
information about particular operator /state. 

When the solution is found the energy can be computed from 

where e a (z) is a single magnon energy 5 It is a simple multi- valued function and ph/mir 
indicate its branches (see Appendix A for notations). Similarly Y^ h and Y^ w correspond to 
the different branches. Finally Zj are exact Bethe roots satisfying the exact Bethe equation 

Y^( Zi ) = -1 . (1.5) 

In this paper we will focus on the sl(2) subsector of the correspondence. From the gauge 
theory point of view this subsector is a closed subsector, with operators being composed 
from scalar fields Z and covariant derivatives T>, 

Tr (V M Z L ) + ... . (1.6) 

The dots stand for all possible permutations of the derivatives with Z fields. In the Y- 
system language, M is the number of roots Zj whereas L enters through large a asymptotic 
of Y# a . From the string theory side, this subsector describes strings moving in AdS^ x S 1 
contained in the larger AdS§ x S 5 space-time. At infinite coupling the string motion is 
classical while the first 1/VA corrections are the quasi-classical one- loop effects which we 
study in this paper. For this subsector the solution is s — ► — s symmetric Y a<s = Y a> ^ s , 
which makes our consideration more transparent. 

In the limit when the number of operators L goes to infinity, Y^ 1T are exponentially 
suppressed and the integral term in ( |1.4|) becomes irrelevant. The remaining term is the 
sum of the individual energies of magnons. In this limit the solution of Y-system is known 
for arbitrary coupling g = ^ and arbitrary state Q. The finite size effects enter in two 



different ways: firstly, the second term in (L4) becomes important and, secondly, the Bethe 



roots Zi move away from their asymptotic values giving rise to the modification of the first 
term. 



4 We strongly believe that Hirota equation associated with this Y-system should automatically imply 
(1.3). From this point of view Hirota equation seems to be more fundamental. See also discussion in |J. 

Magnons are the fundamental spin chain excitations around the BPS vacuum tr Z L . They are dual to 
the string worldsheet excitations. Their dispersion relation and the scattering matrix describing magnon 
scattering are known at any coupling and are properties of the infinite volume theory. 
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In this paper we solve the Y-system equations ( |1.1| ) in the strong coupling limit. We 
construct explicitly the Y-functions for large number of Bethe roots M ~ g 1. Then we 
compare the results with the quasi-classically quantized string using methods developed in 



23, 24, 25, ^(j, 27, 2S, 30 1 and show their perfect match. We show that the Y-system resolves 
the know disagreement of asymptotic Bethe ansatz with the semi-classically quantized 
strings due to finite-size effects. 

2. Y-system in the scaling limit 

The scaling limit is the strong coupling limit g — > oo where the string can be described 
classically. In this limit, the number of Bethe roots M and the operator length L go to 
infinity as g. We also assume that the Bethe roots Zi ~ 1 and \zi\ > 1. This limit of Bethe 
equations was introduced in [£| . The Bethe roots form continuous cuts in the complex 
plane z. They become the branch-cuts on the classical algebraic curve which we describe 
in Sec. |3.1| . 

In the scaling limit we can neglect i/(4g) shifts in the arguments in l.h.s. of (|l.l|) which 
with 1/g 2 precision becomes a set of algebraic equations. This is the key simplification 
which will allow us to completely solve the Y-system at strong coupling. In the notations 
( |1.2j ) we have three infinite series of equations 



Yol = (l + Yo^Xl + Yo^J , s = 3,4,..., (2.1) 

(i + i/Y Aa+1 )(i + i/Y Aa _ 1 ; 

(1 + Y A J 2 
(l + l/Y. a+1 )(l + l/Y -b _ 1 ) 



^ = (1+1/^X1+; - - a --' M — (2 ' 2) 

Yl = T, , J 1+ wf J2 ^ x , « = 2,3,..., (2.3) 



plus four equations 

2 (1 + Ye)(l + Y % ) 

A2 (i + i/Ya 3 )(i + i/y 8 ) ' l ' j 

2 _ (l + Y Qj )(l + Y) 

^ " (1 + 1/ Y) ' {2 - b) 

2 _ (1 + Y Q! )(1 + Y % ) 

^ " (1 + 1/Y A2 ) ' (2 ' 6) 

Y 2 - £ + ^ 2 (2 7) 

^ -(! + !/«,)• (2 ' 7) 

From the asymptotic solution of which we re-consider in the next section, one 
expects that Y-functions have several branch cuts going from z = ±1 + for some 
integers n, to infinity. The approximation we use above is only accurate far enough from 
the branch-cuts. Close to the branch-cuts even a small shift in the argument could cause 
a large jump of Y's. This means we can safely use (2.1-2~7t|) above and below the real axis 



and on the interval [—1, 1] of the real axis, but not on the whole real axis. In this section 
we will mainly focus on the spectral parameters with —1 < Rez < 1 for which (2.1-2?/]) 
are valid. 
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Equation ( |l,3|) also simplifies in the scaling limit 

fl(+)f?(-) = il x(z) - xj l/x(z) - xj i 
R (-) B (+) 11 x( z ) - x+ l/x(z) - ~ /(*)/(«) ' 

where 



= z + iy 7 ! - z 2 , »j = x ph (zj) = Zj + ^zj - + 1 (2.9) 

and xj = x ph (£j ± j-) for the standard definition of y / TT7 with branch cut along negative 
part of the real axis. The functions f(z) and f(z) are defined in terms of the resolvent 



1^1 x 2 



^ X — Xj xj — 1 

in the following way 

/(*)=exp(-iG(x(*))) , /» = exp( + iG(l/x(z))) . (2.11) 

The kernels TZ^ 0a ^ and ^^ 0ct \ defined in Appendix A simplify dramatically in the strong 
coupling scaling limit. For example lZ^ a \z, w) — B^ 0a \z,w) ~ 5(z — w) so that equation 
fll.3| ) becomes 

1 oo 

F = Y®Y@ = T7- II + y O • ( 2 - 12 ) 

0=1 

One should keep in mind that the possible distribution of the Bethe roots is constrained 

M x + M 

V-log^-V I 3 =27rm + 0(l/g) , m€Z, (2.13) 



by 



which reflects the cyclicity symmetry of the single trace operators. 

We see that in the scaling limit we have to solve an infinite set of algebraic equations 
(2.1-2.7, 2.12 ). In the next section we construct the asymptotic solution of these equations 
for Y% a <C 1 and expand it in the scaling limit. 

2.1 Infinite length solution at strong coupling 

In this section we will study the Y-system in the asymptotic limit l# a <C 1 and then expand 
it at strong coupling. We will see that in this scaling limit the asymptotic solutions of Q 
can be recast in the following simple form: 

* w - (s -Af-i, «. _ (sr ; +1 : ^ _ _ , {2 . 14) 

where 

1 f f(f-l) 2 f 



- 5 - 



with ( 2.11 ). The middle node Y- functions are equally simple: 



Ym 



•a 



/(/ - irr - (f - iyff a 
//(/-/) 



where 



A = exp 



Tg + 27TmZ 



(2.16) 



(2.17) 



governs their exponential suppression for large L/g. We now derive these equations starting 
from the asymptotic solution constructed in H in terms of the eigenvalues of SU(2\2) 
transfer matrices for representations with rectangular Young tableaux T as . We have 



l + VY 0s {z) 



T ljS (z + i g )T hs (z 



4s ' 



l + Y A (z) 



T a Az + i-)T a>1 (z 



— ) 



. (2.18) 



-i(z)T 0tS+1 (z) ' ' ' ' T a _ 1; i(z)T a+1)1 (z) 

The middle node functions Y% a (z) are suppressed but can be also expressed in terms of the 
transfer matrices using 



Y. a {z) = {T a , x ) 



*(*- 3j) 



r a— 1 



4g / 



11 11 



; 26-l 

' 4 9 



(2.19) 

where fiW, flW are defined in the Appendix A and a(z,Xj) for Imz > coincides with 
the Beisert-Eden-Staudacher dressing phase |1(| and is understood as an analytical contin- 
uation otherwise. The SU (2\2) transfer matrices for symmetric (Ti jS ) and antisymmetric 
(T 0j i) representations can be found from the expansion of the generating functional [34, 35] 
(see |36|] for some mathematical details). For the s[(2) subsector it reads 

-2 



BM(z + ^)RM(z-^) , 
5<-)( 2 + £)fl(-)(*-£) 



!-(-!)' 



4S ' 



1.9 ' 



1 - D 



(2.20) 



where D = e ' 2 s isa shift operator. One should expand ( 2.20| ) in powers of D and commute 
them to the right 



oo 
s=0 



49 



D s 



W- 1 = J2(-l) a T t 



al Z 



a=0 



;a-l 



D a . 



(2.21) 



We want to construct Y as using this construction in the scaling limit. At strong coupling 
the shift operator D serves just as a formal expansion parameter D, since the shifts it 
creates are suppressed. In the notations of the previous section W becomes 



[l-Df/f)(l-D) 



[1-D/fY 

Now it is very easy to find the general expression for T\ s and T a \ 

,af(f-l) 2 f a -f(f-l) 2 f a 



( a -l)//- S (/ + /) + ( S + l) 
■Lis — /.„ j J- al 



(-1) C 



//(/ " f)f a 



(2.22) 



(2.23) 
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which implies the above mentioned expressions ( [2.14 ) for the Y- functions. In the expression 



for the middle node Y's (pl9|) it is enough to use the leading strong coupling expression 
for the dressing phase - the Arutyunov-Prolov-Staudacher (AFS) phase || 



where x^ = x(z ± j-). Expanding this expression at strong coupling for — 1 < z < 1 we 
obtain Y % ~ (/ 2 A) a T^ 1 , or ( ^16|) . 



We end this section with an important comment which will be later used in Sec. 2.5 



We notice that for \z\ > 1 the second factor in ( |2.19| ) is exponentially small even for finite 
L/g (notice that x(z) denotes the mirror branch x(z) = z + i\J\ — z 2 ) 

I x(z - 




<z + f Q ) " x(z) 



aL 



(2.25) 



and thus the asymptotic solution is accurate for \z\ > 1 even if for — 1 < z < 1 it is signifi- 
cantly modified by the finite size effects. Let us stress that the exponential suppression we 
are discussing is much stronger than the usual finite size exponential suppression at strong 
coupling. The latter is suppressed for large L/g whereas ( p. 25 ) is suppressed for large L 
even if L/g is finite and small. 

In the next section we will analyze the large a and s limit of these asymptotic y's and 
argue that the same asymptotics should be used even when the fine size effects are strong. 

2.2 Boundary conditions 

In this section we propose the boundary conditions which should be used to make the solu- 
tion of the Y-system unique. For that we study the asymptotic large L solution considered 
in the previous section at large a or s and argue that the exact solution should have exactly 
the same behavior. 



From ( 2.14 ) we see that Y^ a oscillates with a because T(z) = f(z)/f(z) is a pure phase. 
To have a well defined large a limit we shift the argument by — iO then \T(z — i0)\ > 1 and 
we get 

a^oo a J 

Similarly 

lim logY. a (,-.0) =lo A/2 

Whereas Y/\ a and l# a decrease exponentially with a, Yq s behave as s 2 . The general solution 
of (|2.1[) with polynomial asymptotics is 

Y 0a = {s-A(z)f -1 (2.28) 

for some A(z). When A is small and l» a are suppressed A(z) is given by its asymptotic 
value ( |2.15| ), otherwise it is some unknown function. In the next section we find its exact 
expression as a function of A. 
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Conditions flgjgp » (1^ ) and (gTgg ) could be seen to be consistent with the TBA 
equations for excited states proposed in |1| (see also |H]]) Y^ should satisfy 6 



log Y Aa (z - iO) = A4 am *log(l + Y^ )-K a _i, m _i*log(l + Y Am )-K a -i®log(l + Y )+log 



P (+) R (+) 



where A4 and X are some kernels defined in Appendix A. At strong coupling the last term 
gives 



1 K> + R + f 

hm - log , w ~ log - 

a-^oo a pi J ) J 



(2.29) 

assuming that the other terms are not growing with a linearly this leads precisely to ( 2.26| ) . 
Similarly, ( 2.27| ) and (|2.28|) could be justified from the TBA equations of [19|. 

2.3 Y-system in T-hook 

In this section we solve Q2.lH2.7p together with ( |2~lp . O ne can achieve a considerable 
simplification of this problem by transforming ( |2.1[]2.7| ) into the Hirota equation. For that 
we rewrite Y as in terms of T as 



Y„ 



T as should satisfy the Hirota equations 



^a+l.a-^o— 1,8 



T a + 1 ,sT a —i ,s + T as +iT L 



a,s—l 



(2.30) 



(2.31) 



from which all ( \2. lH2.7| ) follow. The indices of these T as functions belong to the T-shaped 
lattice (see Fig. 2). Another equivalent representation, which follows from ( |2 . 3 1| ) is 



1 + Y„, 



T 2 

1 n 



-a+l,5-^a— l,s 



(2.32) 



It is important to notice that the choice of T as is not unique for given Y as , there is a 
"gauge" freedom 

Ta,a(z) -> gi (z)[g 2 {z)) a T a , s (z) , (2.33) 

which leaves Y'b unchanged for two arbitrary functions gi(z) and gi(z). 

Below we will see how the general solution of ( |2.31 ) can be constructed for the infinite 
vertical strip, which is the upper part of the T-hook. Then we constrain it by the large 
a asymptotic ( |2.26| ) and ( p. 27 ), and match with Yd s given by ( 2.2£| ) (shown as dark gray 
circles on Fig. 2). As a result all Y's are constructed explicitly in Sec. |2.3.2| for finite L/g. 

2.3.1 Solution of Hirota equation in the vertical strip 

We notice that inside the vertical strip of Fig. 2 the Hirota equation Q2.31 ) coincides with 
the recurrent relation for the characters of SU(4) for the representations with rectangular 



J the — iO is due to the prescription to go below singularities in the convolutions [ p^[ . 



- 8 - 




Figure 2: T-shaped "fat hook". Small black filled circles correspond to the finite T as function. 
Big circles show where solution for y-functions for the strips is applicable. 

Young tableaux. It is well known that the characters are given by Schur polynomials for 
Young tableau (Ai, . . . , A4): 7 

s,{y) = ^f^' 1 -'"- 4 • (2.34) 

de %i )l<M<4 

We only need symmetric solutions T a _ s = T a ^ s which implies that 2/4 = 1/2/1 and 2/3 = 1/2/2- 
For rectangular representation Aj< s+ 2 = 0, Aj >s+ 2 = a, from (|2.34 ) we get T a ^ = 1 and 



This solution of ( |2.31| ) has only 2 parameters because it is suitable for a finite semi-infinite 
strip boundary conditions, appropriate for characters. On the other hand, the most general 
s — > — s symmetric solution should have 4 free parameters. It is clear that we can get one 
more by shifting a in ( [2.35 ) by an arbitrary function. What is also true, although not 



as trivial, is that we can shift independently a in y\ and in 2/2 s ° that the most general 
solution of (|2.31| ) is 



T a ,2 = 1 , 

Ta, = , (^-r + S iyt ) - * ( * + fttf ) ) (2.36) 

(yi - 2/2) (2/12/2 - 1) V v{ - 1 \SiVi ) V2~ l \ S 2V2 

Ti rri'2 rjn rj~\ 

a,0 — J a ,l ~~ i a+l,l- f a-l,l > 

all the others T's are given by T a _ s = T QjS . 



7 We would like to thank V.Kazakov for discussing this point. 
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To establish a relation with the previous section we reparameterize our solution in 
terms of 4 new parameters in the following way 



2/1 



Then 



iVTe- a / 2 



2/2 rj-\ 



S-2 



TU 
~S~ 



T - SUe a+2 ST a+3 - Ue a+2 /T c 



(2.37) 



(2.38) 



'SU{l-T){T -e) V 1-e T 2 — e 

which for small e coincides with ( |2.23| ) up to a gauge transformation ( [2.33 ) and thus zero 
e limit is the asymptotic limit we consider in the previous section. We can now easily 
compute the large a limit of Y 9a and Y& a : 



lim 

a— >oo 



\ogY Aa (z-iO) 



1 logY. a (z-iO) 
log — , hm = loge . 

1 a— >oo a 



From this, using (HH) and ( p7| ), we fix two of the four unknown function, 

= A/ 2 . 



rj-i f 



f 



(2.39) 



(2.40) 



So far we have constructed all l# a and Y\ a in terms of two yet unknown functions S(z) and 
U(z) in such way that ( |2.2| ) and ( p.3j ) are already satisfied. We have to glue this solution 
with the solution inside the right and left wings ( |2.2£| ) parameterized by a third unknown 
function A(z) and demand the remaining equations (2.4f 2.7 , 2.12| ) to be satisfied. 

2.3.2 Matching wings 

In the previous section we construct the solution for the upper wing of y-system (the light 
gray dots on Fig.2) and fl2.28Q gives the general solution for the left and right wings Yd s . We 
parameterized these functions in terms of three unknown functions S(z), U(z), A(z). We 
still have to match these two solutions in the different domains and find the two remaining 
fermionic ^-functions Y®{z) and Y${z). To fix these five functions we have exactly five 
remaining equations ( 2^ - p.7| , ^.12 ). Excluding 1® and J® we get 



y A 2 2 (i + i/yA 3 ) 

1 + V^A 2 



1 + Ym, 



1 

A 2 



(A- 



1)2 
l) 2 



1 



y 2 (i + i/y^) 



F 

(A-l) 2 (F-l) 



where F = YgY®- F can be expressed in terms of Y^ from ( [2,12|) : 

00 

F = J7 + 

J J m=l 



(2.41) 
(2.42) 
(2.43) 

(2.44) 



Notice that the right hand side of (|2.41 - iQ4| ) depends on A and F only while the left hand 
side of these equations does not depend on the two unknown functions S, U and the known 
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functions /, / and A. These equations are relatively easy to solve perturbatively in A. For 
example for A with C(A 4 ) precision we found 

A — j^j + j—^ + A(/ + / - 2) + A 2 (p + f 2 - 2) + A 3 (/ 3 + /3 _ 2 ) + ... . 

We see that the expansion coefficients are very simple. We can easily sum them up to get 
the exact result 

1 //-I 1 9 
A = = h w r - tz h 



A/-1 (/-!)(/-!) /A-l A-l' 
_ 2(/-l) 2 (/-l) 2 (/-l)(5//- 3/- 3/ + !)(/- 1) (/+/-2//) 2 - 
' " /(A-l) 3 / /(A-l) 2 / /(A-l)/ 77 ' 

/(/-1) 2 (/A-1) 2 

/(/-1) 2 (/A-1) 2 ' (A45j 
u _ (/-1) 2 (/-1) 2 



/ 2 / 2 (/A-l) 2 (/A-l) 2 • 

The result is quite simple compared to what one may expect to get from the high degree 
polynomial equations and apparently there should exist some more straightforward way to 
get this result. It is easy to check that fl2.41|,|2.42D and ( ^43|) are indeed satisfied. To check 
( 2.44 ) one can use ( |2.32| ) to get rid of the infinite product 

OO -.00 rp2 -1 rri rp rp 

f =^i na+^o = Ti n T r = ji¥ nm 7^ = ^ A • ( 2 - 46 ) 

// £± x ff ^ ?m-l,oTm+l,0 // To j0 a^oo T a+lfi T 0fi 



which allows to express the r.h.s. as a rational function of e, U, S and T so that (2.45) can 
be easily checked. Notice that we found all Y- functions except for the fermionic ones for 
which we only presented explicitly the form of their product F. However, using e.g. ( |2.7j ) 
we can easily get Y® in terms of the other Y's that we just fixed. 

We can now plug the functions we just found into the Y ayS functions to get explicit 
expressions for all of these functions in terms of /, / and A alone! We recall that these 
functions are completely fixed in terms of the Bethe roots ( |2,ll , 2~l7| ). Since the results are 
not particularly simple, we present them in Appendix B in Mathematica form. 

2.4 Energy and Momentum 

Having all y's computed we can easily evaluate the energy of the state, corresponding to 
a given distribution of roots from Ql.4| ). Using it at strong coupling 

<t(*k) = 4 ± ± + o(X) , er(z) = -^tL=, (2.47) 
x%-l \g z J VI - z 2 

and applying the trick from the previous section to compute infinite products of 1 + Y as 
used in ( |2.46| ) to get 



e"° = fT(l +*.)'= hm (r ^ o)a+1 = (/A-l) 4 (/A-l) 4 

™ T 0i0 (T a+1 ,o) a (A-1)4(//A-1) 2 (/ 2 A-1)(/ 2 A-1) 



(2.48) 
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From where we get the following stunningly simple expression, accurate to all orders in 
wrapping 



Using the expressions from Appendix B one can see that each separate term in (|2.4q ) 
is significantly more complicated then the resulting product! One can easily check ( 2.4§| ) 
using explicit expressions for Y% a from Appendix C and expanding both sides of the equality 
in powers of A. 



In ( 2.4SQ , the integration goes over z G (—1,1) because outside these region Y* a are 



strongly suppressed (see (|2.25|) ). Notice that the first term is of order M ~ g ~ y/\ and 
contains the classical string energy, whereas the second is ~ 1 and should be a part of 
the one-loop correction. The reader may already suspect that the numerator of ( 2.4§| ) 



corresponds to 4 + 4 fermionic fluctuation modes whereas the terms in the denominator 
correspond to 4 - modes of S 5 and 2 + 1 + 1 - modes of AdS§. We make this relation more 
precise in Sec.^. To separate the classical string energy from the one-loop corrections in 
the first term of fl2.52|) we should find the equation determining positions of z% with 1/g 
precision. In the next section we will consider the exact Bethe ansatz equation ( |1.5|) in the 
strong coupling limit. 

The total momentum of the state can be computed similarly to the energy. One should 
simply replace the expression for the magnon energy e a in (|1.4j) by the magnon momentum 



X; 1 f 1 dz 1 



F = ETflT + f ^^j=d z M . (2.50) 
~^ 9{xf ~ 1) 2g J_! 2vr Vl - z 2 



The natural extension of the cyclicity condition ( 2,13| ) is 



P = 2vrm . (2.51) 

This is an additional constrain and one should prove its consistency with the other equa- 
tions. We will assume ( [2.51 ) to be satisfied. 



2.5 Exact Bethe ansatz equations 

In the previous section we saw that the exact energy of a given state in the semi-classical 
limit is given by 

p_V*^+l . f 1 dz z 9 , (/A-l) 4 (/A-l) 4 

U A - 1 l-i 2vr Z g (A - 1)*(//A - 1) 2 (/ 2 A - 1)(/*A - 1) ' {2 - b2) 



and the functions /, /, A determined uniquely in terms of the Bethe roots z% ( |2.11 , 2~T7| ). 



We still need to find the positions of the roots Zi in order to get the exact energy of the 
state. In this section we shall derive that these Bethe equations, accurate to all orders in 
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wrapping with one-loop precision, read 
l M - + 



x exp 
where 



1 = -M: [I— -- 77 _ (*fc,Zj) 

V**/ fJixl-Xj l-l/(x k xj) 

(r(x k , z)M + - r(l/x k , z)M- + u(x k , z)M^jdz - iP 



(2.53) 



M+ _ (/A - l) 2 (/A - l) 2 

" (/2A-1)(//A-1) ' " CPA-1)(//A-1) ' l - j 

The factor a 2 (zk,Zj) contains both the leading Arutyunov-Frolov-Staudacher (AFS) || 
and the sub-leading Hernandez-Lopez (HL) phase |39|. Initially the AFS phase was de- 
signed to give an agreement with classical theory. Then it was realized that an extra phase 
is needed in order to get an agreement with the semi-classical one-loop string energies |58| . 
Basing on the known expressions for the one-loop energies of particular classical solutions 
1 43] this extra phase was found in 3£j. However, in pll] it was shown that even with these 



both dressing factors the Bethe ansatz equations are misses some exponential corrections. 
In [25| l the one-loop compatibility of the asymptotical Bethe ansatz was proven for a 
generic classical string motion in AdS^ xS 5 . In |24], |25| it was also noticed that in order 
to get the agreement one should drop some definite exponential in L/g terms. In the next 
section we will work with finite L/g keeping all the previously dropped terms and show 
that the above equation, obtained from the Y-system, describes accurately the one-loop 
string energies for the s[(2) sub-sector. 

We now derive the above mentioned equations. To find the exact positions of the 
Bethe roots zi one has to evaluate Y 9l (zi) on the physical real axis (see Fig. 3). The results 
we obtain in the previous sections are applicable in the domains where Y's are smooth 
functions, however close to the Bethe roots one may expect poles and the approximation 
used so far is no longer valid. 

To get round these difficulties we use the representation of the y-system obtained in 



I IS], based on the TBA approach for the ground state |20], 19 1. The equation we need is 
the integral equation for the middle node 

i og y % =T lm *iog(i+y % j+2^ 10 )®iog(i+%)+2^ 10 )®^ m _ 1 *iog(i+yA m )+^ , (2.55) 

where Ti m is a kernel containing the dressing phase. * stands for convolution with inte- 
gration over real axis, whereas ® is a convolution with integration around the cut (—1,1) 
see Appendix A for more details. To use this equation we will need to know lA m and 
on the whole real axis. In Sec|2.1| we noticed that for \z\ > 1 the asymptotic solution from 



Sec.2^1 can be used since Y^(z) are strongly suppressed for these values of z whereas for 
\z\ < 1 the solution of the y-system was built above. 

We denote by Y® s the asymptotic solution, constructed for a given set of exact Bethe 
roots Zi in Sec.^^, satisfying the exact Bethe equation Y% l {zi) = —1. The asymptotic 
solution should satisfy the following integral equation 

log y° = 2n( 10 k$ io g (i + y °) + 2iz^m m -i * fog(i + y A °j + i$ , (2.56) 
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Ym a are supressed 



Figure 3: Branch-cut structure of Yu 1 . Dashed line is a "mirror" real axis. Between the branch 
cuts all l» a are exponentially suppressed at strong coupling and one can use the asymptotic solution 
here. To find the exact positions of the Bcthe roots one should analytically continue Y 9l under the 
upper cut to reach a "physical" real axis. 



we can subtract it from ( 2.55 ) to get 



log 
1* 



log ^ = T lm * log(l + Y.J + 27^ 10 )©log + 2K^m m -i * log 

As a result of these trick we do not need any more to go outside — 1 < z < 1 region in the 
convolutions, since the integrands vanish there. We can now analytically continue Y^/YQ 
to the physical real axis where the Bethe roots are situated (see Fig. 3). A similar analytical 
continuation was already performed in || 

|£ = ^ * Ml + *J + 27^°)^® log (2.57) 
+ 2K^ h '™m m -i * log ([±|jr) + 2K m ^(z k - i- g ) * log (i±f^) ■ 

Now we simply have to expand the kernels at large g and substitute Vs. To expand T^ ,mir 
again can use the AFS dressing phase ( p. 24 ) 

U^ h ' mir (z k ,w) ~ r(a*,«;) , fi( 10 ) ph ' mir (z fc ,u;) ~r(l/a*,uO , 

Tim mn i z k,w) ~ -m[2r(x k ,w) + 2u(x k ,w) +p(w)] , (2.58) 
■Km(*fc - w ) — K w - z k )+m [r(x k ,w) + r(l/x k ,w)] , 

where we use the following notations 

x 2 d z 1 x d z 1 , , d z 1 

r(x,Zj = — = -— , U(X,Z) = ttt^ T , P(z) = - — ; , . 

v ' 1 x 2 - 1 2irg x - x(z) ' v ' ; x 2 - 1 2irg x 2 (z) - 1 ' i y 1 4irig ^1 - z 2 
After that we can rearrange the terms in ( 2.57] ) according to kernels r. Using the following 
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Figure 4: Density as a function of z for a typical distribution of bethe roots z% on the real axis. 



"magic" products 8 

l+Y ( l+Y A \ m 1 (/ 2 A-1)(//A-1) 



e +M - 



1 + Vtf fr ( l + YKA - 2 = (/A -I) 2 



i + V% ii 2 V i + J (f 2 A - l) (//A - l 



we get the corrected ABA equations ( |2.53| ) accurate to all orders in wrapping with one-loop 
precision. Here we use the expression for the total momentum ( 2.50| ). Notice that the last 
term in the exponent is irrelevant due to ( p. 51 ). 



2.6 Finite gap solutions 

In this section we expand the corrected Bethe ansatz equation ( |2.53| ), obtained in the 



previous section for a particular type of configurations of roots. Before expanding ( 2.53j ) 



one should take log of both sides. Due to the 2ni ambiguity of the log function one should 
assign an integer mode number for each root x^, i.e. for each of the M equations we can 



assign log(l) = 2iririk in the left hand side of (2.53). For the finite gap solutions we assume 



that the set of mod numbers contains only a finite number of different integers. In this 
limit the Bethe roots xj~ are forming dense distributions along some cuts C n on the real 
axis with some density p{x) so that the number of cuts is equal to the number of distinct 
mod numbers (see Fig.[|). One can establish a one-to-one correspondence between such 
configurations and finite gap classical strings motions considered in the next section. 



Strong coupling expansion of ABA for these configuration was studied intensively [37 



38, 23, p9j. Here we can use the existing results to expand our corrected by wrappings 



Bethe equation ( 2.53| ). Usually one defines the quasi- momenta through resolvent ( 2,1C| ) 



P5(s)= 9 2 7 +G(x). (2.61) 
x z — 1 



8 To compute these products we again use z—iO prescription to ensure their convergence. This prescription 
is inherited from the TBA equation where the integration should go slightly below the real axis. 



From the definition of the resolvent ( |2,10| ) we see that p(x) has poles at x = Xf.. When the 
number of roots goes to infinity these poles condense into branch cuts and we can rewrite 
the quasi-momenta in terms of the density p of the roots 

wW -fe^gj> + W M-**. (2.62) 
x 2 - 1 ^ J Cn x - y y 2 - 1 v ' 

Equation ( |2.53|) gives an integral equation on the density of the roots 

2jm = [p^(x + iO) + p%(x — iO)] + a(x)p'§(x) cotpg(x) + V(x) 

2i I (r(x,z)M+-r(l/x,z)M- + u(x,z)Mojdz , x 6 C n , (2.63) 



where a(x) = g ^-i) ■ ^ ne secon d term in (|2.63| ) is so-called "anomaly" term |37], |38|, |4l|, 23 1 
and the third therm is the contribution of subleading term in the dressing phase - the 
Hernandez-Lopez phase [39| 

v W . a(x) ± 1 '-')'•- \ (%-%), Q. = i * °& . (2 . 64) 

^-^ 7r (s — r){r + s — 2) \ x s x r J J 2m x s 

r,s = 2 
r + s S Odd 

The last line in ( |2.63D incorporates the finite size effects. Finally using the standard nota- 
tions 

±x + G(0) 

Pi (x) = -PzO-lx) , p~ 2 {x) = x2 _ - (2.65) 



we can rewrite Ws in terms of the quasi-momenta 
o M _ 



e 



(l _ e -»P2-*P2) 4 (i - e -m-m) 4 



(1 _ e -2ip^4( 1 _ g-»Pi-jp2)2(l _ e ~ 2ip 2)(l - e~ 2ip i) 

(1 _ e -*P 2 -*P2) 2 

.v, _ (l_ e -^i-^)2 



= ^V4t ' : — T , (2-66) 

(1 _ e - 2 *J>2)(l - e ~ tp i~ ip i) 



(l_ e -a»pi)(l_ e -ipi-*a) 



In the next section we will see how these structures appear in the quasi-classical string 
quantization. 

We have found a set of equations which are supposed to correct the Beisert-Staudacher 
asymptotic equations with the Beisert-Eden-Staudacher dressing phase in the strong cou- 
pling scaling limit. The latter are known to describe the semi-classical string spectrum up 
to exponentially suppressed finite size corrections as described in the previous section. This 
extra corrections which we just derived ought to cure the known mismatch and correctly 
incorporate all wrapping corrections to 1-loop precision. In the next section we show that 
this turns out to be precisely the case! 
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3. Quasi-classical string quantization 



In this section we review the quasi-classical quantization method. Then we consider a 
generic solution of string equations of motion inside AdS% x 5 1 and compute its one-loop 
energy. 

The one-loop correction to a classical string energy could be understood as zero point 
oscillations of fluctuations around the classical solution. To compute it one can expand 
the classical action up to the quadratic order around the classical solution and then find 
the spectrum of oscillation modes. These modes could be labeled by the mode number n, 
which tells us how many wavelengthes fit into the string, and polarization. There are 8 + 8 
bosonic and fermionic polarizations which we label by double indices (ij): 



Bosonic: (1,3), (1,4), (2,3), (2,4), (1,3), (1,4), (2,3), (2,4), 
Fermionic: (1,3), (1,4), (2,3), (2,4), (1,3), (1,4), (2,3), (2,4). (3.1) 

We denote the energies of the vibrations Qn ■ Then the one- loop correction is simply a sum 
of halves of these fluctuation energies (4(1 [hi 



\ J2 , (3-2) 



n,(ij) 

where F{j is +1 for bosonic polarizations and —1 for fermionic. 

The direct computation of these 0,n is only possible in the simplest cases ji^, 43]. 



For a generic solution it is enormously hard to perform this calculation starting from the 
classical action. The tool which allows to handle the quasi-classical string quantization 
efficiently is the algebraic curve technique developed in 23, 26, 27, p9f| . Below we describe 
the construction of the algebraic curve and the method of the quasi-classical calculations. 

3.1 Classical algebraic curve 



The classical equations of motion of the Metsaev-Tseytlin superstring action [^4[ can be 
summarized in a compact form as the flatness condition Q 

dA - A A A = . (3.3) 

for a connection A(a, r; x) which is a local functional of fields depending on an arbitrary 
complex number called the spectral parameter x and taking its values in psu(2,2|4). The 
fact that the classical equations of motion could be packed into the flatness condition is an 
indication that the model is classically integrable. Indeed, we can define the monodromy 
matrix 

M(x) = Pexp * A{x) (3.4) 

A 

where 7 is a loop wrapping the worldsheet cylinder once. The flatness of the connection 
ensures path independence of the spectral data of the super (4 + 4) x (4 + 4) matrix 
M{x). In particular, the displacement of the whole loop in time direction amounts to a 
similarity transformation and we conclude that the eigenvalues of the monodromy matrix 
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are conserved with time quantities depending on the spectral parameter x. We denote the 
eigenvalues of M(x) as 

{e ip i,e ip 2,e ip 3,e ip i\e ip i,e ip i,e ip *,e ip i} , (3.5) 

where p%{x) and p%{x) are so-called quasi- momenta J7|. The quasi-momenta contain 
information about all conserved charges of the theory, in particular the global symmetry 
charges, including the energy E. The eigenvalues are the roots of the characteristic polyno- 
mial and thus they define an 8-sheet Riemann surface. In general these sheets are connected 
by several branch-cuts. The branch points on this surface are the values of the spectral 
parameter x where two eigenvalues coincide and M(x) cannot be diagonalized completely. 
Different classical solutions correspond to different algebraic curves. For many calculation 
the explicit construction of the classical solution in terms of the initial fields entering into 
the Lagrangian is not needed and can be replaced by the corresponding algebraic curve. 
For example the energy can be computed as a simple contour integral 



E 



It is always possible to define the quasi-momenta so that they vanish at x — > oo. Then, 
however, the quasi-momenta should jump by a multiple of 2tt when passing through a cut 

Pi(x + iO) +pi(x-iO) pj(x + iO) +pj{x -iO) 

J - ^ = 2-nn , x £ C . (3.7) 

The quasi-momenta are restricted by the properties of the monodromy matrix M(x). Due 
to super-tracelessness 

Pi +P2 +P3 +Pi =Pl +P~2 +^3 +Pl > ( 3 - 8 ) 

and as a consequence of the special properties of M(x) under x — ► 1/x transformation one 
has 

pj( s ) = - Pi (l/x) , pi(x) = -2-Km-p^l/x) (3.9) 

Pg( x ) = -P3( l / X ) 1 Pl( x ) = +2vrm -pz(l/x) , 

where m is an integer winding number. 

There are also infinitely many points where two eigenvalues coincide but, nevertheless, 
the matrix M{x) can be diagonalized. The two quasi-momenta pi and pj corresponding 
to the coincident eigenvalues have no singularity and differs by 27rn. One can perturb the 
curve by opening a small cut connecting the intersecting sheets of the surface at these 
points. We label these points by an integer n and a couple of indices (ij) 

Pl (xH) - Pj (xH) = 27rn. (3.10) 

One of the nice features of the algebraic curve is the simplicity of visualization of the 
action variables of this classical integrable theory. They are the contour integrals around 
the branch cuts 

2m Jr V x 2 



1 g ) Pi{x)dx . (3-11) 
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In the standard quasi-classical quantization procedure one should assume them to be inte- 
gers. 



3.1.1 Algebraic curve for sl(2) subsector 

The algebraic curve for the string in AdS^ x S 1 was constructed in [45|. In the general 
framework reviewed in the previous section this sector corresponds to the cuts connecting 
P2 withpg outside the unit circle centered at the origin. Automatically, due to the x — ► 1/x 
symmetry (|3.9| ), we will have reflected cuts connecting with p^ inside the unit circle. 



One can easily build the spectral representation for the quasi- momenta [45 



p~ 2 {x) 



where G(x) = J c \_ x dy. We see that these quasi-momenta are exactly those of Sec-fEB] 



-p % {x) 



-Pi(l/x) 
Pl(x) = - 



= p i (l/x) 

Pli x ) 



2S' 



¥ 9 X + G (°) 

x 2 - 1 
+ G(0) 



+ G(x) 



1 



(3.12) 
(3.13) 



The action variables ( |3.11 ) count the number of Bethe roots constituting the cut. In this 
way one established the map between classical solutions and the finite gap configurations 
of Bethe roots [g5| . 

The equation ( |3.7| ) for the st(2) subsector becomes 

p${x + iO) +p${x - iO) = 2im , x<EC, (3.14) 

which is now an integral equation for the discontinuity p(x). 

3.2 Quasi-classical corrections from algebraic curve 

(a) 

Using the algebraic curve it is also possible to find the spectrum of the fluctuations 
around an arbitrary classical solution using the techniques developed in [^3|, 24, 25. 26, 27, 
|29[ |. The perturbations of the given classical solution are reflected in the algebraic curve 
by extra cuts. The small cuts could only appear in the special points of the curve given 
by ( 3.10[ ). The perturbed quasi-momenta differ from the non-perturbed ones by a small 
amount 5pi(x). The minimal size of the cut is restricted in the quasi-classically quantized 



theory by the condition that the contour integral around this new cut ( 3.11 ) is integer. 

From far away the branch points almost merge and the cut looks like a pole with a 
tiny residue 



5pi(x) 



a(x) 



a(x) 



(3.15) 



x - xa 9(x 2 - 1) 

such that ( 3.1 1| ) counts a single quantum. We see that for given n and (ij) the perturbation 
of the quasi-momenta is pretty much restricted and one can compute the energy shift due 
to this fluctuation. This gives precisely In all details this technique is described in 
(see also @). 



1 26, 27 



£ln has two contributions different by their nature. Firstly, the extra small cut carries 
its own energy as we can see from ( |3.6| ) 

2 



n 



>j 



w(a#) , uj(x) 



x- 



1 



(3.16) 
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n=+oo 



-<g o- — -o — -o- ■ ^- ±1 -j TDOO-o — o o 5>- 



Figure 5: Complex plane of the spectral parameter x. Positions of the fluctuations x 1 ^ are shown 
by the gray dots. 

secondly, it deforms others cuts, changing thus their contribution. This second contribution 
we study in the next section. Let us now see the effect of ( |3.16 ) on the one-loop shift (|3.2|). 
We have to compute the following sum 

5Ei loop = D- 1 )^ «') > ( 3 - 17 ) 

n (ij) 

where xli should be found from (|3.10| ). We rewrite this sum over n as an integral 

*El loop = £(-1)^ / 7? cot(7rn) W , (3.18) 

(ij) J 1 

where the contour encircles the real axis. Next, for each polarization (ij) we change the 
integration variable from n to x via ( 3,1C| ) . The integration over n maps to contours which 



encircle the fluctuation positions xli located outside the unit circle U. Then we can deform 
this contour in the x plane to get an integral over the unit circle, centered at the origin 
(see Fig.||) 9 

^-,. = E(-^/ u |^c„t(^)„ W = / ^ W ^„, (3.19) 
W) 

where U + is the upper half of the unit circle and 

g^o = JJ (i _ e -iPi+i P ^ Fij . (3.20) 
(ij) 

We use that 

J2(-lfii(p i -p j ) = 0. (3.21) 

(ij) 



For each (ij) the sum in n in (3.17) is divergent and one should carefully treat the large n's. The one- 
loop shift written as an integral ( J3.19 ) corresponds to a particular prescription of the large n regularization 
which is analyzed in details in |ESj. Whereas from the algebraic curve point of view this prescription is 
absolutely reasonable from the worldsheet action point of view the complete prove is still missing. We 
should notice however that this question is not a particularity of the finite size effects, it is rather related 
to the question of the validity of the ABA at the one-loop level. 
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The product in (g!20|) goes over all 8 + 8 polarizations listed in ( |3.1[) . Notice that from 



( 3.12 ) pi —pj ~ L/g and the integral ( 3.20 ) is exponentially suppressed for L/g large. This 
kind of terms are not captured by the ABA and as a result the ABA can be only used 
when L/g is sufficiently large. 

To our deep satisfaction we notice that for the sl(2) subsector A/q = — A'fo from ( 2.66| )! 



Moreover by changing the integration variable to z = ^ (x + 1/x) we map the integration 
contour to [—1, 1] segment of the real axis and ( |3.19| ) matches precisely with the second 



term of the expression for the energy obtained from Y-system ( 2.52 )! In the next section we 



show how the corrected Bethe equation ( 2.63 ) arises from the quasi-classical quantization. 
3.3 Back-reaction 

So far only the direct contribution of the virtual sea of the fluctuations was computed. 
We have to take into account the back-reaction - the deformation of the quasi-momenta 



close to the cuts of the initial non-perturbed classical curve. In [25| such deformations 
were considered dropping exponentially suppressed finite size corrections. This allowed for 
a precise derivation of the HL correction to the AFS asymptotic Bethe equations. Here we 
will keep all exponentially suppressed terms since we want to derive a set of exact integral 
equations. We split p^ into the part containing all the small virtual cuts and the 
smooth part p~ r . To write down V§ one should take into account x — ► 1/x symmetry and 
some further analyticity constraints such as poles at x = ±1. The basic rule is that each 
fluctuation Vt/i contributes as a pole at x = xli with the residue a(x) on the corresponding 
sheets i and j and also by a pole at l/x % n due to the constraint ( |3.9| ) (see plj for more 
details) 10 



T a(,)y (-1)^ la(x)T (-l) F - ,, X)T (-!)»' 
2 ? \ * fe(^) 2 -l to 1 ''-* 
The second sum goes over all fluctuations starting at p^ 

Bosonic : (2,3), (2,4), Fermionic : (2,3), (2,4), 



(3.22) 



and in the last term, corresponding to the reflected poles, the sum goes over all fluctuations 
starting at . 

nnlrl lisp if'v I \\\ to find t.hp Hisnnntimiitv nf i 



Now we should use ( |3.14| ) to find the discontinuity of p)? r 



pf(x + iO)+pf(x-iO) + 2V i (x)=2Trn, y£C. (3.23) 

We can convert the sum over n in ( |3.22| ) into the integral over x - precisely like we did 
with the energy and then deform the contour to the unit circle 

/ dy ( oc(x) cUVo , sCUv+ lA , . d v Af- \ a(x) . 1,,. . 

^2 = f KJ I \ +a(x)^—t- a (l x)-f \+-U.p' cot ps + -V(x . 

Ju+ 2iri \ x y z — 1 x — y 1/x — yJ 2 1 2 

(3.24) 



3 compare with equation 20 in 
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There is one important difference that now there is an extra pole at x = y caught when 
deforming the contour to the unit circle given rise to the second term (see p4}| ). We denote 



n 

(2i) 



1 



-iPi+iPjyij 



IK' 

(ij) 



and 



V(x) = d> £ 



u+ 27r 



a(x) a(l/x) 



x 



dyijPl +V2-PI -P2) ■ 



(3.25) 



(3.26) 



y l/x-y, 

taking into account that d y {p^ + p 2 — pi — p 2 ) = d y (G{y) — G(l/y)) one can see that 
this coincides precisely with the contribution of the Hernandez-Lopez phase [25] in the 
Y-system analysis ( p. 64 ). For st(2) we again have N± = — Ai± and after change of the 
integration variable to z = \{x + 1/x) we get precisely the equation obtained in the Y- 
system framework ( p. 63 ). Thus we established the match of these two completely different 
approaches at the level of equations. 



4. Summary and future directions 

In this paper we studied the finite size effects at strong coupling for strings in AdS^ x S\. 
We attacked the problem from two directions - from the quasi-classical string quantization 



using the algebraic curve techniques j26], ^7], [29], |3C§ and from the recently conjectured 
Y-system ||. We found the same result in both cases thus providing a very nontrivial test 
of the latter. We also derived the corrected expression for the energy of ( 2.52j ) 



M 



E 



E 



xf + 1 



X 2 

1 x l 



1 dz 
1 2^7f 



,d z Mo , e 



Mo 



(/A-l) 4 (/A-l) 



(A-1)4(//A-1)2(/2 A -1)(/2A-1) ; 

(4.1) 

where / and / ( 2.1 1] ) are some simple functions of the Bethe roots x% and A ( 2.17] ) is the 
exponential wrapping parameter. The last integral term is responsible for the finite size 
effects and vanishes in the large volume limit. We also found that the Bethe roots should 
satisfy the corrected Bethe equation 



L M 



a n 



xjl 



l/(x+ Xj 



3=1 Xk 



Xj 1 



q-cr (z k ,Zj) 



x exp 



l/{x k X] 

[r{x k , z)M+ - r(l/x k ,z)M- + u(x k , z)M^jdz 



where and A4- contain all exponential wrapping corrections 



Ma 



(/A - 1) 



(/A - 1) : 



(4.2) 



(4.3) 



CPA - 1)(//A - 1) ' (pA - 1)(//A - 1) • 

There are many interesting directions which would be worth exploring: 

• It would be interesting to make a more direct analysis by solving our corrected equa- 
tions (|4.2| ) for some simple configuration of roots and comparing the solution with 
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the sum of fluctuation energies for the corresponding classical solution. For example, 
it would be very nice to repeat the analysis of 31] for circular strings in AdS^ using 
the corrected Bethe equations. 



It would be also interesting to compare the corrected Bethe equations ( |4.2| ) with the 
conjectured generalized Liischer formula [47] in the strong coupling scaling limit. 



In this paper we focused on strings moving in AdS 3 x S\. From the Y-system point of 
view this is an important simplification because the excited states integral equations 
are only available for this sector [O]. On the other hand, from the string semi-classics 



point of view, following [24, 25 1, the derivation of the corrected Bethe equations would 
be a straightforward task. It would be very interesting to perform this generalization 
and to use it as a guiding principle to construct the Y-system integral equations for 
any excited state. 

• It would also be very important to consider analytically some states which cannot be 
treated bythe scaling limit - like Konishi state. 

• The Y-system conjectured in Q for Aharony-Bergman-Jafferis-Maldacena [^] theory 
recently was supported from the 4-loop perturbation theory |5l| at weak coupling. It 
would be also interesting to make some strong coupling test of this conjecture. 

• Related to finite size corrections but at weak coupling one should reproduce the 
results of @] from the Y-system set of equations. 

• Finally one can try to generalize the approach used here to solve the Y-system at 
finite coupling by bringing it to a couple of integral equations like in |l6| (see [^] for 
some first steps). 

In short, there are many interesting open problems to address related to the exact com- 
putation of the AdS/CFT planar spectrum and many simplifications are to be expected. 
We are getting closer and closer to finding the exact solution to a four dimensional super- 
conformal gauge theory for the very first time. The methods developed here could be also 
useful for a wide range of integrable theories. The quasi-classical quantization probes the 
theories at finite volume and provides an important information about hidden structures, 
such as Y-systems. 
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A. Notations 



There are two distinct possibilities to define x(z) which is solution to x + 1/x = 2z 



x ph {z) = z + y/I^ly/z~+l , x mh (z) =z + iy/l - z 2 . (A.l) 

By default we always choose x = x m ". These two functions coincide above the real axis 
and have the following properties under complex conjugation 

^ph = x P h , x = 1/x . (A.2) 

We also use the notation for the Bethe roots 

x 3 = X v\z,) , xf= X v\ Zj ±^). (A.3) 

The single magnon energy and momentum are 

e a (z)=a+ y - \ , n a (z) = -]og- £- , A.4 

x{z + f g ) x(z-%) i x(z~%) 

depending on which x(z) we are using it could be denoted £a h (z) or e™ T (z). 
The kernels we are using in the integral equations are defined as follows 11 

Agn 



K n {z) 



7r(n 2 + 16g 2 z 2 ) ' 



n— 1 m— 1 

2 2 



Knm(z) = £) E K2 { z + ih ir) > 



n— 1 m— 1 

2 2 



2g ) 



KH w \ = __HL i g ^ ^ ^ ^ + _!£_ log 

2tt2 x(z -f g )- x(w + ^) 4vrz s x ( w - ^) 

Binm) w) ^d^ lo V*(* + %)-*(">-%) + 

{Z,W> 2ni ° g l/x(z-^)-x(w + ^) 4vri ° g x(u»-f^) 



(A.6) 



There are two types of convolutions * and @>. The first corresponds to the usual integration 
along whole real axis whereas the second one is a convolution a long a path going from — 1 
to 1 and then back on another sheet e.g. 

7^ n %log(l + Y 9 ) = J dz log(l + Y 9 ) - £ (n0) log(l + 1/Y® 



L We use rescaled kernels compare to K. ncv, (z,w) = 2gK, M (2gz,2gw). 
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where 1/1© is the analytical continuation of Y® across the cut u £ (— oo,— 1) U (l,+oo). 
The kernel Ti m (z,w) is defined in the following way. For Imz > || and Imw > ^ we 
define it using the usual Beisert-Eden-Staudacher dressing factor a [10 1: 

%mM = loga (x(z + f g ),x(z - f g ),x(w + ^),x(w - % )) 

n-l 

+ (B {lm \z + %^)-rt lm \z + %, W )) . (A.7) 



n-l 



*- 2 



The function Ti m (z,w) has the branch-points at z = ±1 + ^ and w = ±1 + One 
should analytically continue between them in z and w. Defined in this way function has 
four branch-cuts going to infinity in z variable starting at z = ±1 ± ^ and four branch-cuts 
going to infinity in w variable starting at w = ±1 ± ^ [52, 19 1. 

In ( |2.57| ) we also use the notation 7^.( 1 o)ph,mir an( ^ 7^' mir , which means that one should 
take lZ( 10 \z, w) (or Ti m (z, w)) and then analytically continue it in the first argument along 
a path going around the branch point z = 1 + For vS- w \z,w) it simply results in the 
replacement x(z ± ^) — > x ph (z ± ^). 

In the main text we also use the following generalized Baxter polynomials 



They are complex conjugates one of another ^(z) = ifT=(a 



B. Explicit expressions for F's 

One can copy the expressions below directly to Mathematica from .pdf. Here we denote 
d = A, fb = /, Ym[a_] = 3^, Yp[a_] = Y^, Yb[s_] = Y 0s 

sbf={e->f~2 d, T->f/fb, 

S->((f-l)~2 fb(d fb-l)~2)/(f (fb-l)~2(d f-l)"2), 
U->((f-l)~2(fb-l)~2)/(f~2 fb~2(d f-l)~2(d fb-l)~2), 
A->l/(l-d f)+l/(l-d fb)+2/(d-l)+f/(f-l)+l/(fb-l) , 

F->-((2(f-l)~2(fb-l)~2)/((d-l)~3 f fb))-(f+fb-2 f f b) "2/ ( (d-l)f fb) 
-((f-1) (fb-1) ((5 f-3)fb-3 f+l))/((d-l)"2 f fb)-f fb+2}; 

Ym[a_]=(U S(T-1)~2 T e'a(e-T) ~2(T~a(e-T~2) (U S e~ (a+2)-T)-(e-l)U e~(a+2)+ 

(e-l)S T~(2a) T"3)"2)/((-(2 U T e~(a+l)(T-S T~a) (S T~a T-l)) 

+T~2(S T~(2a) (T-l)~2-U T e~a(S T~a-1)~2)+U e~(a+2)(U S(T-1)~2 e~a 

-T(S T~a-1)~2)) (-(U T e~(a+4)(S T~a T~2-l)~2)-U T~3 e~(a+2)(S T~a T~2-l)~2 

+2 U T~2 e~(a+3)(S T~a T-l) (S T~a T~3-1)+IT2 S(T-1)~2 e~(2 a+6) 

+S T~ (2a) (T-l) "2 T~6))/.sbf; 
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Yp[a_]=-(((e-l) (e-T~2) T~a(-(U T e~(a+3)(S T~ (a+1) -1) ~2) -U T~3 e~(a+l) 
(S T~(a+l)-l)~2+2 U T~2 e~(a+2)(S T~a-1) (S T~(a+2)-l)+tT2 S(T-1)~2 
e~(2 a+4)+S(T-l)~2 T~(2 a+4) ) ) / ( ( (e-T~2) T~(a+l)(U S e~ (a+3) -T) - (e-1) U 
e~(a+3) +(e-l) S T~ (2 a+5))(U e~(a+2)(T-S T~a)+e T(S T~(a+1)-1) 
(U e~a-T~a)+T~(a+2) (S T~a-T))))/.sbf ; 

Yb[s_]=(s-A)"2-l/.sbf ; 

Y22=(d~2 f4(fb-l)~2(d fb-l)~2- d f ~3 (d(d+l) f b~2-4 d fb+d+1) (fb(d(2 fb-3) 
-l)+2)+f~2(d(d(d(d+4)+l)fb~4-8 d(2 d+l)fb~3+2(d(ll-(d-7)d)+l)fb-2 
-8(2 d+l)fb+d+4)+l)-f (d(d+l)fb~2-4 d fb+d+1) (fb(d(2 fb-3)-l)+2)+(fb-l)~2 
(d fb-l)~2)/((d-l)f fb(d fb~2((d-3)d f~2-2(d-2)d f+d+2 f-3)-fb(d(2((d-2) 
d-l)f ~2-3(d-l)d f+2 d+7 f-4)+f )+(d-3)d f~2- 2(d-2)d f+d+2 f+2 fb-3)); 

Yll=F/Y22/.sbf ; 
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